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We have studied the appearance of vortex arrays in a rotating 4 He nanodroplet at zero temperature 
within density functional theory. Our results are compared with those for classical rotating fluid 
drops used to analyze the shape and vorticity in recent experiments [L.F. Gomez et ah, Science 
345, 906 (2014)], where vortices have been directly seen in superfluid droplets for the first time. In 
agreement with the experiments, we have found that the shape of the droplet changes from pseudo¬ 
spheroid, oblate-like for a small number of vortices to a peculiar “wheel-like” shape, delimited by 
nearly flat upper and lower surfaces, when the number of vortices is large. Also in agreement 
with the experiments, we have found that the droplet remains stable well above the stability limit 
predicted by classical theories. 
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Helium-4 droplets created by expanding a cold helium 
gasi or fragmentation of a criogenic liquid attain a lim¬ 
iting temperature below 0.4 and constitute the only 
self-bound superfluid systems. Superfluidity in helium 
droplets was established through the dissipationless ro¬ 
tation of an OCS molecule inside them, as indicated by 
the appearance of a clean ro-vibrational spectrum^ More 
recently, the indirect evidence of quantum vortices^ - — and 
the frictionless displacement of swift impurities in helium 
droplets^ point towards a superfluid character of helium 
nanodroplets. 

Superfluid 4 He droplets cannot be set into rotation as 
ordinary droplets or rigid bodies. If a rotating helium 
droplet in the normal phase above the superfluid transi¬ 
tion temperature T\ = 2.17 K is cooled down reaching 
the superfluid phase, it reacts by storing its angular mo¬ 
mentum either into quantized vortices or into travelling 
capillary waves £ Conversely, a critical angular velocity 
uj c has to be supplied to the superfluid droplet for the 
nucleation of vortices with quantized velocity circulation 
in units of h/M , where h is the Planck constant and M 
is the mass of a 4 He atom. Single vortices in helium 
droplets have been addressed theoretically by methods 
of different complexity, see e.g. Refs. I9M131 

When the angular velocity is increased above cj c , larger 
amounts of angular momentum may be stored into the 
superfluid by increasing the number of nucleated vortices. 
These vortices arrange themselves into ordered structures 
(lattices) whose existence in bulk superfluid 4 He was es¬ 
tablished long agO i 14 i 15 We refer the reader to Refs. 16- 
lll for a general presentation of the subject. 

Very recently, superfluid He nanoscopic droplets in 
fast rotation have been studied by coherent X-ray 
scattering^ The existence of vortex lattices inside the 
droplets was established by the appearance of Bragg pat¬ 
terns from Xe clusters trapped in the vortex cores in 
droplets made of N = 10 8 —10 11 atoms (corresponding to 


radii from 100 to 1000 nm) produced by the fragmenta¬ 
tion of liquid helium expanding into vacuum. The shapes 
of the droplets were consistent with those of axially sym¬ 
metric oblate pseudo-spheroids with large aspect ratio 
(AR), defined as the ratio of the long half-axis length b 
to the short half-axis length a along the rotational axis. 
While normal liquid drops change their shape as rota¬ 
tion becomes faster— to resemble a “peanut” (multi¬ 
lobe shape) or a “blood cell”, no evidence of such shape 
shifting has been seen in the helium nano droplets*^ As 
shown in the following, this is fully confirmed by our cal¬ 
culations. 

The presence of dopants was instrumental for detect¬ 
ing the vortex cores, although their number was sensi¬ 
bly smaller than the number of helium atoms ( Nx e ~ 
10 -3 7VHe) and their presence is not expected to intro¬ 
duce large deformations in the droplet in spite that they 
locally distort the superfluid around them, see e.g. Refs. 
H0ll23l . Possible effects on the distribution of vortex cores 
inside the droplet might come from the additional rota¬ 
tional energy associated to the Xe mass, especially at the 
periphery of the droplet. Although such effects seem to 
have been observed occasionally in the experimental im¬ 
ages of Ref. E we will not consider them here. In the 
case of a rotating nanocylinder, these distortions were 
found to be negligible^ 

Once experimentally established the presence of a vor¬ 
tex lattice in a droplet of aspect ratio b/a , the number 
N v of vortices in the lattice could be determined ap¬ 
proximately from the vortex areal density (Feynman’s 
formula^) 

n v = N v /S = 2Muj/h (1) 

where uj is the rotational angular velocity and S = tt b 2 is 
the equatorial cross section of the droplet. Since uj can¬ 
not be directly determined in the experiment^ the anal- 
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ysis relies on the classical relationship between the AR - 
experimentally accessible through the diffraction contour 
maps- and the angular velocity, whose connection with 
the parameters a and b is given by the classical theories 
of rotating liquid drops. 20,21 In this way the experiments 
estimated that the number of vortices in a single droplet 
could be as large as N v = 160. When the vortex density 
is particularly large the experimental images showed also 
the occurrence of “wheel-shaped” droplets^ which have 
no classical counterpart. 

The distinct features of superfluid helium, namely its 
irrotational flow and the possible appearance of quan¬ 
tized vortices, are of course not included in the classi¬ 
cal rotating droplet model . 2Q i 21 The existence of a large 
vortex lattice might influence the appearance of the ro¬ 
tating droplet, and the irrotational moment of inertia is 
known to be very different from that of the rigid body^ 
These facts call for theoretically addressing rotating he¬ 
lium droplets with accurate methods which have proven 
to provide reliable results for superfluid 4 He in confined 
geometries. 

We present here a Density Functional Theory (DFT) 
study at zero temperature of pure superfluid helium 
droplets hosting an increasing number of vortices. To 
our knowledge, this is the first realistic study of multi¬ 
vortex configurations in 4 He nanodroplets. A previous 
attempt to study multi-vortex configurations in super¬ 
fluid droplets is described in Ref. |26|, where a simpli¬ 
fied model assuming linear vortices and a rigid spherical 
droplet was used. 

We have recently analyzed a simpler model system, 
namely a rotating superfluid 4 He nanocylinder hosting 
arrays of linear vortex linesthat constitutes the start¬ 
ing point of the present study. Within our approach, a 
self-bound superfluid 4 He droplet is described by a com¬ 
plex effective wave function T(r, t) related to its atomic 
density as p(r, t) = |T(r, t)| 2 . In the fixed-droplet frame 
of reference (corotating frame) we seek for stationary so¬ 
lutions T(r, t) = e - ^/^<F(r), where the chemical poten¬ 
tial /I and the time-independent effective helium wave 
function <F are obtained by solving the time-independent 
equation 

[H-wL z ]${ r)=/i$(r), (2) 

where H is the DFT Hamiltonian, 27 L z is the angular 
momentum operator around the z- axis, and u is the an¬ 
gular velocity of the corotating frame. 

To determine <F(r) describing a configuration where 
N v vortex lines are present we follow the “imprinting” 
strategy, i.e. we start the imaginary-time evolution of 
Eq. 0 leading to the minimum energy configuration 
with a helium wave function^ 



FIG. 1: Helium droplet configurations hosting (from bottom 
to top) N v = 2,4, 7, and 9 vortices. The left column shows 
the density in the z = 0 symmetry plane (top view), while 
the right column shows side views (x = 0 plane). 


where po(r) is the density of the vortex-free droplet and 
(, Xj,yj ) is the initial position of the j- vortex linear core 
with respect to the z -axis of the droplet. During the 
functional minimization the vortex positions will change 
to provide, at convergence, the lowest energy vortex con¬ 
figuration. It is worth stressing that we work in Cartesian 
coordinates and that no symmetry is imposed to the so¬ 
lutions of Eq. ©. We refer the reader to Ref. [23| and 
references therein for technical details on how this equa¬ 
tion has been solved. 

Due to the high computational cost of our calculations, 
we have limited this study to a helium droplet made of 

_ 1/3 

A^He = 15000 helium atoms having a radius R = roN Ue ' 
with ro = 2.22 A, i.e., R = 54.7 A. This droplet is still 
much smaller than the experimental ones, which in turn 
limits the number of hosted vortices. However, our find¬ 
ings can be compared with the experimental results on 
much larger droplets once scaled with a dimensionless 
characteristic rotational velocity D defined as21 


$o(r) = y/po( r ) 

3 =i 


(x - Xj) + i(y - yj) 
yOr-Xj ) 2 + (y-yj ) 2 


( 3 ) 


n = 


lMp 0 R 3 


87 


( 4 ) 






























3 



FIG. 2: Circulation lines of the velocity field corresponding 
to the N v = 9 configuration of Fig. [T| The inset displays in 
a color scale the regions around the vortex cores where the 
modulus of the velocity held is higher (bright spots). The 
dark spots are regions of low vorticity due to the interference 
between the velocity fields of neighboring vortices. 


where po = 0.0218 A 3 is the helium atom density and 
7 = 0.274 K A “ 2 is the surface tension of the liquid. 
For the Nn e = 15000 droplet, Q = 1 corresponds to 
oj = 1.13 x 10 10 s -1 . 

Figured] shows configurations hosting N v = 2,4, 7, and 
9 vortex arrays obtained with Q = 0.43,0.54,0.62, and 
0.69, respectively. Comparing the top and lateral views, 
it is apparent that the droplet becomes increasingly de¬ 
formed, oblate-like, as N v (and thus £4) increases. Also 
apparent is how the droplet surface locally deforms and 
the vortex lines bend forced by the physical requirement 
that their open ends hit perpendicularly the surface. The 
bending is smaller for larger N v and, at variance with the 
classical droplet results^ - — 2 the droplet becomes “wheel¬ 
like” as indeed observed in the experiments^ 

By increasing the angular velocity the number of 
vortices that can be stabilized inside the droplet in¬ 
creases. Eventually, a maximum number of vortices can 
be hosted, above which the rotating droplet will no longer 
be stable. For the Nn e = 15000 droplet we have found 
that the maximum N v value is 9. 

As in rotating buckets^ the higher the angular ve¬ 
locity, the more packed the vortex array is around the 
rotation axis. This leaves a “strip” around the equator 
of the droplet free of vortices that can be clearly appre¬ 
ciated in the N v = 9 case as shown in Fig. [2] where we 
display several circulation lines of the superfluid velocity 
field. The inset shows in a color scale the regions around 
the vortex cores where the modulus of the velocity field is 
higher (bright spots). As expected, the calculated circu¬ 
lation of the velocity field of the superfluid along a path 
surrounding the vortex array equals N v , and equals unity 



FIG. 3: Stability diagram for a number of vortex lines N v = 
0,1, 2,..., 9 as a function of the dimensionless angular veloc¬ 
ity fh The zero of the energy scale corresponds to the energy 
of the vortex-free droplet. The vertical axis is the energy per 
atom in the corotating frame referred to that of the vortex-free 
droplet. The triangles mark the crossings between different 
stability lines. 


around every single vortex. 

Figure [3] shows the calculated stability diagram. As 
for the rotating bucket ^ 23 1 28 1 29 the energetically favored 
structures for N v > 5 are made of a ring of vortices en¬ 
circling a vortex at the center of the droplet. 

It is worth observing that Eq. ©> which strictly ap¬ 
plies to an extended vortex triangular (Abrikosov) lattice 
made of a large number of vortex lines, is also fulfilled in 
the present case in spite of the limited number of vortices. 
This occurs in the case of N v = 7, where the equilibrium 
structure (see Fig. [TJ), is a “patch” of a triangular lattice 
whose areal density is n v = 2/(a/ 3d 2 ), d being the mean 
inter-vortex distance. By equating this expression to Eq. 
© -with the value Q, = 0.62 used to obtain the 7-vortex 
configuration shown in Fig. [T|- one gets d = 28.3 A. An 
average vortex-vortex distance d = 28.2 A can be esti¬ 
mated from Fig. [lj which compares very well with the 
result of the classical vortex theory^ - — 

Figure [T] shows that, disregarding the vortex array, the 
shape of the droplet is almost axially symmetric. To 
determine its AR we have calculated a and b from the 
moments of the density distribution obtaining^ b/a = 
[(x 2 )/(^ 2 )] 1 / 2 , where ( x 2 ) = fp(r)x 2 dr and (z 2 ) = 
f p( r) z 2 dr. 

The AR dependence on the angular velocity Q is shown 
in Fig. 0J together with the curve derived from the clas¬ 
sical model for a rotating liquid droplet^, and used in 
Ref to interpret their data. 
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FIG. 4: Aspect ratio b/a as a function of the dimensionless 
angular velocity H. The solid line shows the curve obtained 
from the classical model for axisymmetric rotating droplets. 



FIG. 5: Calculated ratio I u /I z shown as a function of the 
number of vortices N v 

The figure shows that for a given angular velocity, the 
classical droplet model overestimates the calculated as¬ 
pect ratio. Most likely, the calculated points in Figs. 0] 
should get closer to the classical curve for larger droplets 


having many vortices, which unfortunately are beyond 
the current possibilities of the DFT approach. Notice 
also that in the experiments of Ref. axially symmetric 
stable droplets were observed with aspect ratios as high 
as b/a = 2.3 corresponding to = 0.71, considerably 
larger than the shape instability threshold of classical 
droplets leading to multi-lobe configurations, Q = 0.56. 
Our calculations also yield a similar behavior. 

It appears from Fig. 0]that as N v increases the depen¬ 
dence of the AR on D within the corresponding stability 
region (i.e. within each group of crosses shown in Fig]4]) 
becomes increasingly important, i.e. the droplet is more 
easily deformed. Such increase of the AR proceeds by the 
flattening of the droplet as the vortex cores are pushed, 
as the frequency is increased, towards the center of the 
droplet. 

Another interesting difference between classical and su¬ 
perfluid behavior, which is likely related to the deviations 
from classical theory just discussed, emerges if we look 
at the ratio between the moment of inertia around the 
z-axis, I z calculated from the droplet mass distribution, 
and that obtained from the response of the superfluid to 
rotation, I u = (L z )/lj. The ratio I u /I z is shown in Fig. 
[5] as a function of N v taking for If a value in the middle of 
each stability region. One may notice that the higher the 
angular velocity the closer the moment of inertia becomes 
to the rigid-body moment of inertia. 

To summarize, within DFT we have shown that the 
shape of rotating helium droplets hosting a number of 
vortices evolves from spheroidal at low angular veloci¬ 
ties to wheel-like at high angular velocities. On the one 
hand, multi-lobe configurations present in classical vis¬ 
cid droplets^ are hindered by the appearance of vortex 
arrays whose regular distribution is hard to accommo¬ 
date into a peanut-like (or higher lobe number) shapes. 
On the other hand, the physical requirement that the 
ends of the vortex lines hit perpendicularly the droplet 
surface favors their parallel alignment for large vortex ar¬ 
rays, and hence the appearance of wheel-like shapes, as 
indeed observed in the experiments. Finally, in spite of 
the apparent differences between normal and superfluid 
rotating droplets, the classical relationship between the 
aspect ratio and the angular frequency is fairly fulfilled, 
the classical relationship underestimating the actual an¬ 
gular frequency by less than 10 % for the relevant, larger 
vortex arrays. Thus, it can be used with some confidence 
in the analysis of the experimental results. 
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